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1 Introduction. 



In this paper we continue our investigation of 7* — 7* scattering at high energies (see Ref.pj 
for our previous attempts to study this process in the DGLAP dynamics). We concentrate our 
efforts here on the case of two photons with large but almost equal virtualities. It has been 
argued |21 E! that this process is the perfect tool to recover the BFKL dynamics jlj which is the 
key problem in our understanding of the low x ( high energy) asymptotic behavior in QCD. 

It is well known that the correct degrees of freedom at high energy are not quarks or gluon 
but colour dipoles IHl El E] which have transverse sizes and the fraction of energy z. 
Therefore, two photon interactions occur in two successive steps. First, each virtual photon 
decays into a colour dipole ( quark - antiquark pair ) with size rt. At large value of photon 
virtualities the probability of such a decay can be calculated in pQCD. The second stage is the 
interaction of colour dipoles with each other. The simple formula ( see for example Ref. jOj ) 
that describes the process of interaction of two photons with virtualities Qi and (^2 ( ~ Qi ) 
is (see Fig. [T] ) 

Nf 

(t{QuQ2,w) = J d\j2 (1-1) 



a, 6 



^ dzi J rfVi,t|^^_^(Qi;^i,ri,j)|2 ^ dz2 J rfV2,t|^'r,L(Q2; ^2, r2,t)|^ 2 Ar(x, ri,^, ra,*; 6*) 

where the indexes a and b specify the flavors of interacting quarks, T and L indicate the 
polarization of the interacting photons where r j denote the transverse separation between quark 
and antiquark in the dipole ( dipole size) and Zi are the energy fractions of the quark in the 
fluctuation of photon i into quark-antiquark pair. N{x, ri j, r2,t; bt) is the imaginary part of the 
dipole - dipole amplitude at x given by 

X = ^1 + ^2 /^2) 

W^ + Qj + QI 

for massless quarks {W is the energy of colliding photons in c.m.f.). bt is the impact parameter 
for dipole-dipole interaction and it is equal the transverse distance between the dipole centers 
of mass. 

The wave functions for virtual photon are known ^lOj and they are given by (for massless 
quarks) 

I^T(Q;^,n)r = Y.^-^zli^' + i^-^f)Q'Kl{Qn), (1.3) 

a " 

|*2(Q;^,n)P = Y.'^-^zlQ-z\i-zfKl{Qn)-. (i.4) 

a ^ 

with Q\ = z{l — z)Q'^ where Za denote the faction of quark charge of flavor a. 

Since the main contribution in Eq. ()1.1|) is concentrated at ri t ~ VQi ^ V/^ 
ri t ~ I/Q2 ^ Va* where /i is the soft mass scale, we can safely use pQCD for calculation of 
the dipole-dipole amplitude in Eq. ()1.1|1 . 



2 



Figure 1: The picture of interac- 
tion of two photons with virtu- 
ahties Qi and Q2 larger than a 
"soft" scale. 

¥(r , Q^) 

In this paper we study this process in the region of high energy and large but more-less 
equal photon virtualities {Q\ ^ Q\ ^ 1//^^) the framework of the BFKL dynamics. 
In the region of very small x (high energies) the saturation of the gluon density is expected 
|llllT^ IT^. We will deal with this phenomenon using Glauber-Mueller formula 00111 which is 
the simplest one that reflects all qualitative features of a more general approach based on non- 
linear evolution [HI |T2l 1121 El • For 7* —7* scattering with large but equal photon virtualities, 
the Glauber-Mueller approach is the only one on the market since the non-linear equation is 
justified only for the case when one of the photon has larger virtuality than the other. 

In the next section we discuss the dipole-dipole interaction in the BFKL approach of pQCD. 
The solution to the BFKL equation, that describes the dipole-dipole interaction in our kinematic 
region, has been found jT^I and our main concern in this section is to find the large impact 
parameter (fej behavior of the solution. As was discussed in Ref. [1^1 El El P, we have 
to introduce non-perturbative corrections in the region of ht larger than 1/2 where is 
the pion mass. We argue in this section that it is sufficient to introduce the non-perturbative 
behavior into the Born approximation to obtain a reasonable solution at large ht. 

Section 3 is devoted to Glauber - Mueller formula in the case of the BFKL emission (31. 
Here, we use the advantage of photon - photon scattering with large photon virtualities, since 
we can calculate the gluon density without uncertainties related to non-perturbative initial 
distributions in hadronic target. We consider the low x behavior of the dipole-dipole cross 
section and show that the large impact parameter behavior, introduced in the Born cross 
section, fulfills the unitarity restrictions ( unitarity bound JHl)- Therefore, we confirm that 
the large ht behavior can be concentrated in the initial condition (see Refs. [121 El Q] without 
changing the kernel of the non-linear equation that governs evolution in the saturation region 
as it is advocated in Ref. [T7j . 

In the last section we summarize our results. 

2 Dipole-dipole interaction in the BFKL approximation. 

In this section we discuss the one parton shower interaction in the BFKL dynamics ( see Fig.E}. 
We start with the Born approximation which is the exchange of two gluons (see Fig. |21 ) or the 
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Figure 2: Dipole-dipole interaction in the Born approximation. 



diagrams of Fig. El without emission of a gluon. 



2.1 Born Approximation: 

These diagrams have been calculated in Ref. [T] using the approach of Ref. [201 and they lead 
to the following expression for the dipole-dipole amplitude: 
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{rx,ur2,tM) = nag ^,,3 In -^^ — ——^ -— - (2.5) 



2^c \ {b - ZiTi - 2:2^2)^ - ZxTx - ;Z2^2)^ 

,2 ^2 



2^c -1 1 2 /PI,1'P2,2' 



where is the fraction of the energy of the dipole carried by quarks; Zi = Zi — \ and pj^^ = 
Pi~pk- pi is the coordinate of quark i (see Fig.|21). AH vectors are two dimensional in Eq. ()2.5|1 . 

Each diagrams in Fig. |2l is easy to calculate j2l] and the first diagram is equal to 



, 2 ^^c 1 2 i„ 2 



Summing all diagrams we obtain Eq. (j2.5|) . 



lnpi,i' lnp2,2' • (2-7) 



N(ri,r2;x,b) = 




One parton shower 




V(ri,Ri;V) 

(O(v)y 

V(r2,R2;-v) 



BFKL Pomeron 



BFKL emission 

Figure 3: One parton shower interaction in the BFKL approach. 



We are interested mostly in the hmit of large bt ^ ri ^ 
amplitude can be reduced to a simple form. 



r2,f where the dipole-dipole 



- 1 rl. 



iV2 

c 



{21 



after integration over azimuthal angles. 

Therefore, we have a power-like decrease of the dipole-dipole amplitude at large ht , namely 
]\[BA ^ "[IjiZ^a, Such behavior cannot be correct since it contradicts the general postulates 



of analyticity and crossing symmetry of the scattering amplitude ^H]- Since the spectrum of 
hadrons has no particles with mass zero, the scattering amplitude should decrease as e"^'"''** 
[T^ . In Ref . P we suggested a procedure of how to cure this problem which is based on the 
results of QCD sum rules |i21j. Following this procedure we rewrite the dipole-dipole amplitude 
as the integral over the mass of two gluons in t-channel; and we assume, as in QCD sum rules, 
that this integral describes all hadronic states on average. Restricting the integral over mass 
by the minimal mass of hadronic states ( 2 ) we obtain the model which provides the 
exponential fall at large ht ^ 1/(2 m^r) and does not change the power like behavior for small 
ht < 1/(2 m^. 

We choose for the Born amplitude the following formula 



iV''^(ri,t,r2,t;6t) = vra; 



N^c 3 

One can easily see that Eq. ()2.9|) reproduces Eq. ()2.8p and leads to 

2 - 1 r\^ m% 



i^4(2mA) • 



iV''^(ri,t,r2,t;6i) 
at large 6t ^ 1/(2 mj^). 



TT 



iV2 



2m^ht 



(2.9) 



:2.10) 



2.2 BFKL equation: 

The emission of a gluon is described by the BFKL equation |3] which was solved in Ref. [T3] for 
fixed ht (see Ref . [221 123 1211 many useful discussion of the different aspects of the solution). 
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The solution can be presented in factorized form (see Fig. |21). 

N{x,r,^t,r2,t;bt) = (2.11) 

J Zm 

with 

= ^ (2^(1) - - + ; (2.12) 
where ip{f) = dlnr{f)/df ,r{f) is Euler gamma function and where 

V(rit,Ri;u) = — "-^ (2.13) 

using the following notations: y = ln(xo/x); r^.f is the size of the colour dipole "z" and Ri is 
the position of the center of mass of this dipole. In Eq. (j2.1ip function (j)in{y] ?"2,t; bt) should be 
found from the initial condition which determines the dipole amplitude at fixed x = Xq, namely, 

A^(x = xo,ri,t,r2,t; 6t) = A^-^^(x = xq, ri,t, r2,t; 6*). 

It should be stressed that the BFKL equation is a linear equation in which the kernel does 
not depend on ht (see Ref. [13]). Therefore, (piniy] h) could be an arbitrary function on bt- 

In Eq. fl2.11|) we can take the integral over R2 which leads to 



■ + 1 1/ 



V(r2t,R2;-i^) = I — ^ ^ 1 • (2.14) 

We are interested in the large bt behavior, namely, bt ^ ri^t ~ '^2,1- It is instructive to 
consider two cases: 

• DLA: 0. This is so called double log approximation of pQCD (DLA) 

in which we consider j -C r2,t and Q;s'ln(l/x) ln(r|j/r^J 1 while ^^^(l/x) ^ 1 
as well as as^'^^ir'^t/r'lt) ^ 1 a,nd as <^ 1- We have considered this case in Ref. pLj and 
found that the emission of gluons does not induce any additional dependence on bt which 
is concentrated only in the Born amplitude. Indeed, we can see this property directly 
from the solution of Eq. ()2.11|) . 

Integrating over Ri we find that the integrand of this integral falls down rapidly for 
Ri > bt due to Ri dependence of the vertex V{r2,t, R2]^) (see Eq. fl2.14|) ) providing a 
good convergence for the integral. For Ri < bt we can neglect Ri dependence of the 

vertex V{r2^t, R2', ^) and consider it as ( j . The integral over Ri of V{ri^t, Ri, i^) 



for Ri < bt gives (ri ^ 
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Therefore, J d'^RiV{ri^t,Ri,i^) V{r2,t,bt,u) (^t/^?) Kt/^i)^"*''- Finally, taking 

4>in{i^', ri t] ht) = Ti as ("^tt T ('^i t ^tY Ki{2 111^, bt) 1 \ . , the dipole amplitude has a 
form 

N^'^^{x,n„r2,t;h) = N^\x,n,,r,,-A) I ^ ^^i'^)y + ih-^'^)HrlJrl.) (2.15) 

J 111 I 

Considering r2,t ^ '"i,* and taking into account that io{v) 1 -iu at | ~ ^ ^ 

one can take the integral in Eq. ()2.15p explicitly. The answer is well known ( see Ref. 
for example), namely, at low x 



iV^^^(x, ri,„ r^x, h) = N^^ix, r,,; h) I, |^2^ ^ V Hrlt/rl,) | (2.16) 

for fixed coupling constant^. 

Diffusion approximation: f <^ 1. For such small values of v the integral over 
Ri is convergent for Ri > ri t (see Ref. ) and, therefore, we neglect the Ri dependence 
in V{r2,t, h, J^)- Introducing a new variable ^ = Ri/ri^t we see that 

Jd'R,V{n„R,;u) = {rl/^-'"' Jd'^(i+^nni-^nry^'^ (2.17) 

where n is a unit vector in the direction of r^^j. The integral is a function of u only and 
can be absorbed in (pmiv] h) in Eq. ()2.11|) . For V{r2^t, h, at bt ^ r2^t we have 



2 \ 5+*'' 



V{r2„h,-u) = (2.18) 



Therefore, the dipole amplitude is 



,2 ^2 \ h+i" 



N^^{x,n,,r2,uk) = J ^^ct>U^-A) e-^^^y ■ (2-19) 

We choose (i)in{i'] h) to be of in the form 



A^^ - 1 . 1 
0„(i/; bt) = nas -f-p^ {ni^ hf K^{2 h) . (2.20) 



At small values of v we can expand ^^{i') 



uj{u) = ul - Du^ (2.21) 



§In this paper we consider only the case of fixed QCD coupling since the BFKL equation is not proven for 
running as- 
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with 

= ^^4 1n2; D = UC{3) ; (2.22) 

TT TT 

Finally, we can evaluate the integral over u in Eq. ()2.19|) using the method of steepest 
decent and obtain the following expression for dipole amplitude: 

N''^{x,n,t,r2,t;bt) = (2.23) 



In 



2 2 
2 ^l,t ^'2.t 



^ ( n,t r2,t ml hi )K,{2m^h) e^^ ^ " 

At 1/(2 m^) > bt > n^t ~ r2,t 

N''^{x,n,t,r2,t;bt) ^ (2.24) 



2(iv,^-i) n,tr2,t nr 



2 2 
, 2 ''l.t ''2.t 
In • — T 



while at 6^ > 1/(2 m^) 

iV''''(x,ri,i,r2,t;6t) ^ (2.25) 



2 2 
2 ^l.t ^2,t 



In^ 

^'^L y TtT7, 2 6t 

3Ar2 \} 2 ni^bt^ Dy 



3 Saturation and unitarization in the Glauber - Mueller 
approach. 

3.1 Glauber - Mueller formula. 

The Glauber - Mueller approach El El takes into account the interaction of many parton 
showers with the target as is shown in Fig. |3] In our case of more or less equal but large 
virtualities of both photons this approach gives a unique opportunity to study the high energy 
asymptotic behavior of the dipole amplitude since other methods based on non-linear evolution 
equation [TTl |T21 UHl Ull 12^1 do not work in the case of two dipoles with more-less equal 
sizes. 

The main idea of this approach is that the colour dipoles are the correct degrees of freedom 
for high energy scattering (this idea was formulated by A.H. Mueller in Ref. |S|). Indeed, the 
change of the value of the dipole size (Ar^) during the passage of the colour dipole through 
the target is proportional to the number of rescatterings (or the size of the target R) multiplied 
by the angle kf/E where E is the energy of the dipole and kf is the transverse momentum of 
the t-channel gluon which is emitted by the fast dipole. 
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Figure 4: The Glauber - Mueller 
approach for the dipole-dipole 
scattering amplitude. 



A r+ oc R — . 

E 

Since kt and rt are conjugate variables and due to the uncertainty principle 

1 



(3.26) 



h oc 



Therefore, 



kf 

Art oc R— < rt if i? < rl E oi x < 
E 



2m R 



(3.27) 



Since the colour dipoles are correct degrees of freedom , they diagonalize the interaction 
matrix at high energy as well as the unitarity constraints, which have the form 



2 N{x,ri^t,r2,t]ht) 



\aei{x,ri,,r2±]ht)? + Gin(x,ri,,r2uk 



(3.28) 



where Qei is the elastic amplitude of the dipole-dipole interaction and is the imaginary part 
of a^i {N = Imaei). 

Assuming that the amplitude is pure imaginary at high energy, one can find a simple solution 
to Eq. (IT^ . namely 



aei{x,ri^t,r2,t;bt) 
Gin{x,ri,,r2,ubt) 



(1 



(3.29) 
(3.30) 



where Q is the arbitrary real function. 

In Glauber - Mueller approach the opacity Q is chosen as f2 = 2 N'-^^^{x,ri^t,r2/, ht) where 
]\fOPS jg ^YiQ dipole-dipole amplitude for one parton shower interaction that has been found in 
the previous section (see Eq. 1)2.111) ). 



3.2 Saturation. 

One can see that if we substitute the explicit solution to the BFKL equation of Eq. ()2.23p at 
any fixed bt the opacity Q = 2N^^ increases at x — > 0. Therefore, the dipole-dipole amplitude 
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given by Glauber-Mueller formula of Eq. (|3.29p tends to unity in the region of low x. This 
statement is called saturation fTR IT^ IT^ since the physical interpretation of is the density 
of colour dipoles at least when N is not very large. In this discussion the saturation appears to 
be the consequence of unitarity for fixed bt. However, we have learned several examples where 
the dipole density could reach a maximum value without having any effect on the elastic dipole- 
dipole amplitude at fixed bt (see Ref. [in] and paper of Kovchegov and Mueller in Ref.jESl)- 
However, for 7* — 7* scattering of two small dipoles the initial condition is given by Born 
amplitude of Eq. ()2.5|) which is small. Therefore, we have no reason to expect that the dipole 
density will be high due to the final state interaction. 

3.3 Unitarization. 

To obtain the unitarity bound for the dipole-dipole cross section we have to integrate over bt, 
namely 

a {dipole - dipole) = 2 j dHt N^^^ {x,n,t,r2,t;bt) = 2 J dHt {l - e-^°^(^'^i.*''-2,*;fe0^ (3.31) 

Following Froissart [TH] , we divide the region of integration over bt in Eq. ()3.3H) in two parts 

a{dipole - dipole) = 27T " db] N^^\...-bt) + / db] N^^\...-bt) (3.32) 

where bQ{x) is defined from the equation 

iV^^(x,ri,t,r2,t;6o(x)) = 1 (3.33) 

It is easy to see that for bt < bo{x) N'^'^'^ < 1 since A^-^-^ > 1, while for bt > 1 and for 
A^-^-^ < 1 A^*^^ < A^-^-^ Therefore, we have the following unitarity bound 

a{dipole - dipole) < 2nlbl{x) + H db^^ N^^\...;bt)] (3.34) 

V Jbo{x) J 

Let us consider two possibilities. The first one that ^0(2;) <^ 1/(2 m^). In this case the 
solution to Eq. ()3.33|) follows directly from Eq. ()2.24p for the amplitude A^^^ and for y 1 

ln(^^) = (3.35) 

or 

blix) oc n^tr2,t e'^y (3.36) 
Substituting Eq. ()3.3(ij) into Eq. ()3.34|1 we can obtain 

a{dipole - dipole) < 2 7Tbl{x){l + 2} oc eV^^^ (3.37) 
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where the second term is calculated by integrating first over bt Eq. (|2.19p and after that using 
saddle point approach. Since ^saddle turns out to be small at low x and we neglect it. 

Therefore, in this kinematic region we face a power -like increase of the dipole-dipole cross 
section as was pointed out in Refs. [TTj. 

However, this power-like increase will stop for 60(2;) > 1/2 m^r. Indeed, for such large values 
of ht we should use Eq. fl2.25|) for the dipole amplitude N^^ . For such large values of bo^x) 
Eq. has a solution which at low x is 

= + 0{lny) (3.38) 



which leads to 



2m,, 



a( dipole — dipole) < 27cbl(x) = — ^ ln^(xo/x) (3.39) 

4m£ 



which comes from the first term in Eq. (j3.34p . It is easy to understand that the second term 
in this equation gives a term which does not increase with y. Eq. (j3.39|) is the unitarity bound 
which has the same energy dependence as for hadron-hadron collisions [19^ but in our approach 
we are able to calculate the coefficient in front of y^. The bound of Eq. ()3.39p is the same as 
was derived in jTHl HE] • 

It should be stressed that the diffusion approximation that we used was derived only at 
small values of saddle point in z/ integration in Eq. (j2.19|) which is equal to 

In ''^f 

WsaddleI = Jf!^ « 1 (3.40) 
2Dy 

at bt = boi^x) from Eq. 



3.4 Saturation scale. 

Eq. ()3.33|) does not have a solution at any values of ri t and r2,t (formally, we obtain a negative 
values of bo{x)). The same equation at 6f = 0, namely 

N{x,r,,t,r2,t;h = 0) = 1 , (3.41) 

determines the saturation scale. At ^ > Vgat the opacity Q in Glauber-Mueller formula is 
larger than unity [Q > 1), Eq. ()3.33p has a solution and we are in the saturation region with 
Eq. (I3.39P for the unitarity bound. If ri t < Vgat, opacity Q < 1 ai any value of bt- This is a 
domain of perturbative QCD in virtual photon scattering. 

N{x,ri^t, r2,t] bt = 0) we can find from Eq. ()2.11|1 integrating over Ri, namely 
N{x, n^t, r,^t; bt = 0) = I ^ 0.„(z/; h = 0) rf^ i?i e'^^^^^ V{n^t, Ri, ^) V{r2,t, Ri, -v) ■ (3.42) 
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Since Vgat from Eq. (|3.4ip is much smaller than r2,t we need to find Eq. (j3.42j) only for ri ^ <C r2^f 
This observation simplifies the calculations. Indeed, the main contribution in the integral 
over Ri stems from Ri <^ r2^f Therefore, we can neglect the Ri - dependence of vertex 
V{f2,t, Ri] which has the form 



16 X^^^'^ 



V{t2,uRv^-v) = y-^j . (3.43) 

To perform the integration over Ri we use the following formula ( see equation 3.198 of Ref.^T] 
) 

B{- - iz/, - - iv) I — ; ^-^ I = (3.44) 

^ ^ . . ,,l- ./_r, , . / 



dt (t(l-t))-^^M i?? + (l-2t)i?i ■ ri,, + ^ 



where B{ji,u) = r(/i)r(z/)/r(/i + u) is the Euler beta - function. Integrating Eq. (|3.42|) over 
Ri using Eq. ()3.44|) we obtain that 



iV(x,n„r,,6, = 0) = /^0a-;&Oe-(^)^4f^-^4^^f^y"^"" (3.45) 



The Born approximation at 6^ = and at r2^t ^ '"i,* can be reduced to PP 

N^\n,, rxu h = 0) ^ nal . (3.46) 

^^c ^2 ^2 '2,t 

It is easy to choose (pmiy] h = 0) in such a way that the final answer for N{x, ri^t, 1^2/, h = 



is: 



Nix, n„ r,, b, = 0) = .al f ^ ^ e^^'"' " "^^-'^-^ • (3-47) 

Z2 Z2 J 2 711 ^ — lu 

We can find the solution to Eq. 1)3.411) in the saddle point approximation for the integral 
over z/ in Eq. ()3.47p fn\ 1^ . Introducing a new variable 7 = | + zi/ we have the following 
equation for the saddle point value of 7 = 75: 



I7=7S 



y+Hrit/rL) =0 (3.48) 
Substituting Eq. in Eq. we obtain 

N{x,ri^t,r2,t;bt = 0) oc g'^^^^^-^s^'^HtKJ = e?/M7s) - 7s ^k^.s) (3.49) 
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Using Eq. (j3.49|) we can solve Eq. (j3.4ip in semiclassical approximation ( see Ref. |2H] in 
which we cannot calculate the numerical factor in front of Eq. ()3.49p . Indeed, N{x, ri r2/, bt = 
0) is constant on the line 

^l.=7o y + HrLix)/rl) = (3.50) 
with 7o is the solution to the equation [HI EB] ' 

7o d-f 



17=70 



(3.51) 



Eq. ()3.5Up leads to a power-like increase of the saturation momentum {Qsat{x) = '2,/rsat ) at 
high energies ( low x). Namely, 

"(70) "(70) 

QlM) oc ^ (-) " - Ql (-) " (3.52) 

rl^t \xj \xj 

Actually, the pre-exponential factors in the steepest decent method of taking integral over 7 
could change the x-dependence of the saturation scale adding some log(l/x) dependence in 
Eq. (|3.52p ( see Ref. j2H] for an analysis of such corrections) . 



3.5 Unitarity bounds for 7* — 7* scattering. 

To obtain the unitarity bounds for 7* — 7* scattering we need to substitute the unitarity bound 
for dipole-dipole cross section (see Eq. ()3.39|) ) into Eq. (jLip and to perform integrations over 
Ti^t and Zi. J cP Ttl"^ l{Q, ZyTt)]"^ is convergent while / ci^ rt|\l/r(Q, z, rt)p has a logarithmic 
divergence that we need to deal with. Eq. ()3.39j) holds only for ri^t > fgat since if ri^t < fgat 
dipole-dipole cross section is small and proportional to / d'^ htN^^^ . As has been mentioned 
we consider ri_t < r2.t- On the other hand Ki{z) ~ 1/z aX z < 1. Finally, one can see 

/ d'n / dz,\^i^T{Q,z,n)\' = CQ-HQlJx)rl,) , (3.53) 

where Cq = Ea ^ ^a- We recall that g2^^( x) r| J does not depend on '"2,t- 

The integration over r2^t is concentrated at the limits ^ < r2^t ^ I/Q2 and leads to 

r 4 1 
/ d^r2,tdz2\'^T\^ = Cq - In . 

^7o was called ko in Ref. 11 and Ao in Ref. 28 . The numerical solution of Eq. H3.51|l leads to 70 — 0.63. 
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Finally, for 7* — 7* cross sections we have 



f * *N ^ ^2 16 / Qsatix) Qsatix)\ ( UJl 2/ / x\ /gcr.N 

-.(7-7*) < (m^) (i^in^(Vx)l ; (3.55) 

9 r gi i l2m2 



-.,.(7* -7*) < ^.Y^ln:^J (^l-^(Vx)J ; (3.56) 
a.,47*-7*) < (^^A^o/x)] ; (3.57) 



Since the saturation scale increases as a power of (1/x) one can see that the energy behavior 
of the unitarity constraints is 



aT,T(7*-7*) < [Cif) r^'lJ'^^' I lAxo/x); (3.58) 



""^^y^ I ln^(^o/x) ; (3.59) 

.2 16\ / TTUJl \ , 2/ 



-,l(7*-7*) < (.C'Qyj (^^jln^(W^)- (3.60) 

Note only (Tl,l{i* — 7*) has the same energy dependence as hadron-hadron collisions ^H] 
but even this cross section has a different coefficient in front. Cq as well the as numerical 
factor 2/3 come from the photon wave function while ul reflects the BFKL dynamics making 
Eq. ()3.60|1 and Eq. ()3.57|) quite different from the unitarity bound for hadronic reactions. 



4 Conclusions. 

In this paper we use 7* —7* scattering as the laboratory for studying the large impact parameter 
behavior of the amplitude in the saturation region. At first sight, this processes occurs at short 
distances for both photons with large virtualities and could be calculated in perturbative QCD. 
We demonstrated that the non-perturbative QCD corrections have to be introduced for large 
bt even for this process. The main result of this paper is the statement that it is enough to 
include the non-perturbative QCD corrections in the Born approximation and neglect them 
in the kernel of the BFKL equation. This result confirms the mechanism suggested in Refs. 
[inillHI but it contradicts the arguments of Ref. jTTj. 

This result does not mean that the BFKL kernel correctly describes the large bt behavior. 
The uncertainties in the large bt tail of the kernel will not affect the high energy asymptotic 
behavior of the dipole amplitude. Let us assume that kernel of the BFKL equation can be 
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written as K + AK where K is normal BFKL kernel in pQCD and AK includes the non- 
perturbative contribution. We know that AK oc e"^'^'"''* from general properties of the strong 
interaction JH]- Let us treat AK as a small correction and calculate the first digram of the 
order of AK (see Fig. EJ. 




Figure 5: The diagrams for the 
first order correction with respect 
to Ai^jWhich includes the non- 
perturbative QCD contribution 
at large values of the impact pa- 
rameter. 



The sum of all diagrams in Fig. El leads to a contribution 

AK (l - N^^^{x,ri^t,r2,t,bt)) = AK e-^'''''^''''''"''''''''^ (4.61) 

Since for bt < is very close to unity, the above corrections are suppressed. Only 

for bt > bo{x) we can expect a considerable contribution. However, this contribution is pro- 
portional to e-'^™'^''o(x) — ^-ujLin{xo/x) ^ Therefore they turn out to be very small. 

This simple discussion shows why the strategy to include the non-perturbative corrections 
in the Born amplitude, works 

Actually, the main result of this paper, namely Eq. ()3.38|1 . is based on a simple physics (see 
Ref. JE])- We have demonstrated here that the multi rescattering processes embraced by the 
Glauber-Mueller formula lead to a different resulting parton cascade than is given by the BFKL 
approach. The principle difference is the fact that the multi parton shower interaction creates 
a new scale or mean parton transverse momentum ( saturation scale) given by Eq. ()3.52|) . 

N{x,ri^t,r2,t,bt) denotes the parton density, consequently the fact that N{x,ri^t,r2,t,bt) — >■ 
1 can be understood as the fact that the partons reach a maximal density at low x. This 
phenomenon is called saturation |1HI12[ IT^. Therefore, at low x we have the parton distribution 
in the transverse plane presented in Fig. [HI the uniform distribution of partons ( dipoles) with 
sizes of the order of l/Qsat{x) in the disc of radius R{x). If one of the dipole inside of the disc 
will emit one extra parton this emitted parton will interact with others partons and as a result 
of this interaction its transverse momentum will be of the order of Qsat{x)- It means that this 
emitted gluon will not change its position in impact parameter space since due to uncertainty 
principle 

AbtPt ^ 1 (4.62) 

its Abt ~ i/Qsat{x). However, for the parton at the edge of the disc the situation is dif- 
ferent since the emitted parton in the direction outside of the disc can move freely without 
any interaction. This parton changes the size of the disc by its displacement in bt, namely 
Abt ~ 1/Pt ~ 1/2 m^r. In this estimate we consider the non-perturbative emission with 
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'.'r-tr^Gsverse j)lane j 



Figure 6: The structure of the parton cascade of the fast photon in the frame where the second 
photon is at rest. The picture is the three-dimensional one since the thickness of the vertical 
line reflects the value of the transverse momentum. The thicker the line the larger value of the 
parton transverse momentum. The left part of the picture shows the parton distribution in the 
transverse plane. 

Pt ~ 2m^ because, as have been discussed, a non-perturbative emission is needed to provide 
the unitarization of our process. Since the emission that leads to a growth of the disc occurs 
in one direction ( the exterior of the disc) it leads to i? = < \Abt\ > n where n is the number 
of emission at given x. Since the emission takes place at the edge of the disc where the parton 
density is rather small, N{x, ri t, r2,t, bt) is determined by the BFKL dynamics only JHinHl- In 
the BFKL approach j^n = cui ln((xo/a;) since N'^^^ oc e" = e*^' iri(xo/x)^ Therefore, we obtain 
Eq. dSSHl), namely, R{x) = b^ix) = ^ ln(xo/x). 

We have discussed in this paper the structure of dipole-dipole interaction in the Glauber- 
Mueller approach which is the only one on the market for the interaction of two dipoles of 
the same sizes. However, for two dipole with small but different sizes the non-linear evolution 
equation fT| ^1 El should be solved to which the BFKL emission is only an approximation 
in the region of small partonic densities. Comparison of the result of this paper with the dipole- 
dipole interaction in , so called, double log approximation shows that the BFKL dynamics 
does not change physics at large bt- The non-linear evolution equation at fixed bf was solved 
PU] in the case when the BFKL kernel was replaced by the double log one. The solution leads 
to the answer in the saturation region with geometrical scale j2Sl EDI CHI 

N{x,r,,t,r2,t;bt) =f{t= r?,, QL(a^) e"""^-''* ) . (4.63) 

Therefore, we believe that for the BFKL dynamics Eq. ()4.63p will hold. This belief is based 
on the similarity between double log and BFKL approximation for 7* — 7* processes. 
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A Appendix 

The integration over Ri in Eq. fj2.11|) can be taken explicitly \2'i\ 1^ and Eq. ()2.11|) can be 
reduced to 

SRiV{n^uRi\y)V{n^u\Ri - ht\]-v) = ,i , ^ (A.l) 



li + 



ci x*^ F{h, h, 2h, x) F{h, h, 2h, x*) + C2X^" V(l-/i, 1-h, 2-2h, x) F{l-h, 1-h, 2-2h, x*) 

where F{a, P,'y, z) is the hypergeometrical function (see Ref. ); x is the complex 

anharmonic ratio: 

{b - Zi n^t - Z2 r2,t) {b - Zi n^t - Z2 r2,t) 

and h=^+ih'. XX* gives 

XX* = (A.3) 

{b - Zi fi^t - Z2 f2,t)2 {b - Zi fi^t - Z2 r2,t)2 



One sees that Eq. ()A.3|) is invariant with respect to rotation in the plane. 

The coefficients ci and C2 have been calculated in Ref. [23 and they are equal: 

c. = rr2 f^T^y (A-4) 



V r(2/^) ; Vi^-h)) 



(A.5) 



However, one can see that Eq. (jA.lll does not reproduce the Born term of Eq. ()2.5|1 at 
y = 0. To understand why it is so we should consider the vertex V{ri^t,Ri] i^) in momentum 
representation ( see Ref. |2S1), namely, 

nri,i,Q;z/) = J d^Rie'^V{n^t,Ri;i^) . (A.6) 
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BFKL Pomeron 






Figure 7: Structure of the Mueller- Tang vertex. 



It turns out 23^ that 



4 



4 



4 
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At small Q Eq. ()A.7|) leads to the following behavior of vertex V{r — 1, t, Q\ v) 



/vi2 / ^2 



As have been discussed the matching with the Born approximation occurs at iu 
limit 



(A.7) 



(Ai 



|. In this 



(A.9) 



which has correct analytical behavior. Actually this behavior dictates the choice of the coeffi- 
cients Ci and C2 in Eq. ()A.4jl and Eq. ()A.5jl . 

However, at —iu | the low Q behavior has a singularity Therefore the symmetry 

of Eq. 1)2.111) with respect to sign of u is broken. Mueller and Tang [21] pointed out that this 
problem can be cured by adding to the expression of V{ri^t, Q] ^) of Eq. ()2.13p . namely, 



F(ri,4,i?i;z/) - 



(A.IO) 



(^1 + h ri,t) 



{Ri - ln,ty 



As was found[15^ such terms can be added due to gauge invariance of QCD. In momentum 
representation (see Eq. ()A.6|) ) V'^^'^{ri^t,Q', i^) can be written as a sum of three terms as it is 
shown in Fig. [7| 
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The Mueller- Tang vertex leads to the Born approximation amplitude in the form of Eq. (j2.5|) . 
However, as was discussed in Refs. |15| 122 12^; it has not been proven that this vertex will 
satisfy the BFKL equation. The solution in the form of Eq. ()A.1|) has a different form of the 
Born amplitude, namely. 



However, these two expressions for the Born amplitude are equivalent due to gauge invariance 



Using Eq. ()A.1|) we can calculate the dipole-dipole amplitude at 6j = and, therefore, the 
saturation scale with better accuracy than within Eq. (j3.47|) . On the other hand the saturation 
momentum increases for x —>■ . Such an increase guarantees that Eq. (j3.47|) approaches the 
amplitude given by Eq. ()A.1|) in the region of low x. This is the reason why we prefer to use a 
simple solution of Eq. ()A.1|) instead of full expression of Eq. flA.l|) . 

It is easy to show that Eq. ()A.1|) describes all properties of diffusion approximation that has 
been discussed in section 2. 




(A.ll) 



of QCD T5\ 
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